Tzitzeica solitons vs. 
relativistic Calogero-Moser 3-body clusters 
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Abstract 

We establish a connection between the hyperbohc relativistic Calogero-Moser sys- 
tems and a class of soliton solutions to the Tzitzeica equation (aka the Dodd- 
Bullough-Zhiber-Shabat-Mikhailov equation). In the 6A^-dimensional phase space 
of the relativistic systems with 2N particles and antiparticles, there exists a 2A''- 
dimensional Poincare-invariant submanifold ilp corresponding to A*" free particles 
and bound particle-antiparticle pairs in their ground state. The Tzitzeica A'^- 
soliton tau-functions under consideration are real-valued, and obtained via the dual 
Lax matrix evaluated in points of fip. This correspondence leads to a picture of 
the soliton as a cluster of two particles and one antiparticle in their lowest internal 
energy state. 
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1 Introduction 

The equation 



6 6 (1.1, 



has a curious history. It first arose a century ago in the work of the Rumanian mathe- 
matician Tzitzeica [1,2]. He arrived at it from the viewpoint of the geometry of surfaces, 
obtaining an associated hnear representation and a Backlund transformation. 



For many decades after Tzitzeica's work, the equation (1.1) was not studied, two 
papers by Jonas [3,4] being a notable exception. Thirty years ago, it was reintroduced 
within the area of sohton theory, independently by Dodd and Bullough [5] and Zhiber and 



Shabat [6], cf. also Mikhailov's paper [7]. In this setting, (1.1) is viewed as an integrable 
relativistic theory for a field ^'(t,!/) in two space-time dimensions, written in terms of 
light cone (characteristic) coordinates, 

t = u — v, y = u + v. (1.2) 



Accordingly, the PDE (1.1) is known under various names, and has been studied from 
several perspectives, including geometry [1-4], classical soliton theory [5-20], and quantum 
soliton theory [21-26]. Moreover, it has shown up within the context of gas dynamics [27, 
28]. 

The principal aim of this paper is to tie in a class of soliton solutions to the Tzitzeica 



equation (1.1) with integrable particle dynamics of relativistic Calogero-Moser type. (A 
survey covering both relativistic and nonrelativistic Calogero-Moser systems can be found 
in [29].) The intimate relation of the latter integrable particle systems to soliton solutions 
of various evolution equations (including the sine-Gordon, Toda lattice, KdV and modified 
KdV equations) was already revealed in the paper in which they were introduced [30], 
and was elaborated on in [31]. Later on, the list of equations whose soliton solutions are 
connected to the relativistic Calogero-Moser systems was considerably enlarged [32-34]. 
In all of these cases, the N solitons correspond to point particles. 

The novelty of the present soliton-particle correspondence is that the Tzitzeica A^- 
soliton solutions at issue correspond to an integrable reduction of the 3A^-body relativistic 
Calogero-Moser dynamics. Physically speaking, a Tzitzeica soliton may be viewed as a 
lowest energy bound state of three Calogero-Moser 'quarks', one of which has negative 
charge, whereas the other two have positive charge. 

A crucial ingredient for establishing the correspondence is the relation between an 
extensive class of 2D Toda solitons and the relativistic Calogero-Moser systems, already 
studied in [33]. Indeed, the relation can be combined with the link between the Tzitzeica 
equation and the 2D Toda equation. The latter link has been known for quite a while, 
and we proceed to sketch it in a form that suits our later requirements. 

Assume that 0„ is a solution to the 2D Toda equation in the form [35] 

(f)n,uv = exp(0„ - 0„_i) - exp(0„+i - (pn), u E Z, (1.3) 
which has the symmetry property 

(p-n = -(Pn, (1-4) 



2 



and which is moreover 3-periodic, i.e., 



Then one has in particular 
so that 



1.5) 
1.6) 

;i-7) 



satisfies (1.1). Conversely, a solution ^ to (1.1) yields a solution </>„ to (1.3) satisfying 
(1.4) and (1.5) when one sets 



0, 



'l.i 



The point is now that there exist soliton solutions to (1.3) that can be made to satisfy 
the extra requirements (1.4)-(1.5), hence yielding soliton solutions to (1.1). The relevant 
2D Toda solitons are those found by the Kyoto school [36,37]. These solitons also formed 
the starting point for [33]. They are most easily expressed in tau-function form, the 
relation of r„ to 0„ being given by 



0„ = ln(r„+i/r„), n G Z. 
In terms of r„, the evolution equation becomes 

dudyhlTn = 1 - Tn-lTn+l/rl, 



and the extra features (1.4) and (1.5) amount to 

T-n+l = 

and 



1.9) 



1.10) 



(1.11) 



Tn+Z = Tn- (112) 

As we show in Section 2, one can make special parameter choices in the 2D Toda 



2A^-soliton solutions r„('U, f) so that they satisfy (1.11)-(1.12). The function 

^ = ln(r2/ri) 



;i.l3) 



then satisfies (1.1), and can be viewed as a Tzitzeica A^-soliton solution. (In Lie algebraic 
terms, the successive requirements (1.11) and (1.12) amount to reductions A^o -Boo 
A 



(2) 



More specifically, the 2D Toda 2A^-solitons of Section 2 are of the form 

Tn{u,v) = det(l2Ar + D{n,u,v)C), 



1.14) 



where the dependence of the 2N x 2N (Cauchy type) matrix C and diagonal matrix D 
on the parameters a,b,^^ G C^^ is suppressed. To satisfy the B^o restriction (1.11) and 
to prepare for the 3-periodicity restriction (1.12), these 6A^ parameters are expressed in 
terms of 2N parameters 0, ^ G and a coupling parameter c. We then show that 
the tau-functions have period / for c equal to ir/l, so that the Tzitzeica restrictions are 
satisfied for c = 7r/3. 
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In Section 3 we make a further parameter change, trading 0i, . . . , for 'positions' 
qii . . . iQn- This reparametrization ensures in particular that the summand involving all 
exponentials has coefficient 1. Restricting attention to the parameter set 

v = {{q,e)eR^'' \eN<---<e^}, (i.is) 

the tau-functions take their simplest and most natural form. In particular, for parameters 
in V the tau-functions are real-valued. Furthermore, their space-time dependence is such 
that the g^'s and 6'j's can be interpreted as relativistic positions and rapidities. Last but 
not least, it is in this form that the c = vr/S Tzitzeica A^-soliton tau-functions can be most 
easily compared to the tau-functions arising in the framework of the 3A^-body relativistic 
Calogero-Moser systems. 

Section 4 is devoted to this comparison. The choice of regime for the Calogero-Moser 
systems is the same as for almost all other soliton equations. Specifically, the regime is 
the hyperbolic one, with the Poincare group generators given by 

N+ N- 

i=i j=i 

S^= exp(±p+)V++ ^M±Pj)Vr, (1.17) 

l<i<Af+ l<j<A^- 

«^)^- n n ;.li-w, ) ■ n-i^ 

l<k<N+,kf!=i ^ 



sinh {xf - xl)/2j \ cosh {xf - Xj )/2 




sin^ c \ -i-r / sin^ c 



sinh2(Xj. -Xi)/2J ^^j^-^^ cosh^(x^. - x+)/2^ 
Here, the generalized positions and momenta vary over the phase space 

n = {{x^,x-,p+,p-) E M2{JV++Ar_) I < ... < < ... < x'}, (1.20) 

equipped with the symplectic form 

UJ = Y dxlAdpl+ Y dxj AdpJ. (1.21) 

l<j<Af+ l<j<N- 



(Only the Landau-Lifshitz solitons of [38] involve the more general elliptic regime 
The coupling c that is needed, however, differs from the value 7r/2 relevant for most soliton 
equations (which correspond to ^i^"*). Indeed, as already mentioned, the connection of the 
space-time dynamics (1.16)-(1.19) to the Tzitzeica tau-functions arises for the '742'-value 

c = 7r/3. (1.22) 

In order to clarify further restrictions and to describe our main result in some detail, 
a few more ingredients need to be introduced. First, the Poisson commuting 'light-cone' 
Hamiltonians 5*+ and S- can be expressed in terms of an (A^4. + A^_) x (A^_|_ + A^_) Lax 
matrix on Q via 

S+ = Tr(L), S_ = Tr{L-^). (1.23) 
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This matrix has a product structure DA4D, with a diagonal matrix D and a matrix 

that arises by suitable substitutions in Cauchy's matrix C{x.y). i.e., a matrix with 
elements (xj — yj)~^, whose determinant is given by Cauchy's identity 

Hence one can explicitly determine the symmetric functions of L and its inverse, and all 
of these functions Poisson commute. 
Specifically, the Lax matrix 

L = DMD (1.25) 



is given by 



D = diag(exp(p+/2)(\/+)^/^...,exp(p+^/2)(\/+Jl/^ 

exp(pr/2)(l^f )^/^ . . . , exp(p^_/2)(l-_)V2), (1.26) 



Mik = 



M 



N++j,N++l 



M 



I sm c 


sinh((a;/" 


-0/2 + 


ic) ' 




i sin c 




sinh((a;~ 


-^r)/2 + 


ic) ' 




sine 





" cosh((x7-4)/2 + ^c)' 
Mi^N++i = -MN++i,i, (1.27) 

where i, /c e {1, . . . , A?".!.} and j, / e {1, . . . , N_}. It is clear from this that L is self-adjoint 
when or vanish. Physically speaking, this is the case in which only particles or 
antiparticles are present. For iV+A^_ > however, L is not self-adjoint. Instead, L is 
J-self-adjoint, that is, the adjoint L* equals JLJ, where 

J = diag(ljv+,-liv_). (1.28) 

When one takes the interparticle distances to infinity, it becomes clear that there is a 
subset of Q on which L is diagonalizable with real eigenvalues. (Of course, for the one- 
charge case this is true on all of Q.) However, already for the simplest non-self-adjoint case 

= iV_ = 1, complex-conjugate eigenvalues are also present, reflecting the existence of 
bound states of a particle and antiparticle. More precisely, choosing from now on 

ce(0,7r/2), (1.29) 

these bound states are encoded by eigenvalues 

exp(^+),exp(^_), §±^e±iK, Ke{0,c]. (1.30) 

Moreover, the case of an eigenvalue pair on the boundary of this sector corresponds to 
phase space points 

x+ — x^ — x, p^—p^—p. (1-31) 
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For p = this yields a 1-parameter set with minimal energy H = 2MoCosc, cf. (1.16)- 



(1.19), and, more generally, no oscillation takes place for the 2-parameter set (1.31) 



The special case just discussed is the only one that can be explicitly understood 
in an elementary way. Already for + = 3 (the simplest case relevant for this 
paper) a complete account involves considerable analysis. The general case has been 
elucidated in great detail in [31], and we need to make extensive use of Sections 2B-6B of 



that paper, which pertain to the c-range (1.29). The action-angle map constructed there 



makes it possible to understand all of the Poisson commuting dynamics at once, including 
their soliton type long-time asymptotics (conservation of momenta and factorized position 
shifts). 

It is beyond our scope to recapitulate even the case + N_ = 3, but we do mention 
the key starting point for the analysis performed in [31]. This is because it clarifies why 
the diagonal matrix 



A{x) = diag(exp(x^), . . . , exp(x^^), — exp(a; 



'1 ;5 • 



-exp(a;^_)). 



(1.32) 



is of pivotal importance. This matrix is referred to as the dual Lax matrix, and its relation 
to L is encoded in the commutation relation 



icot{c){LA - AL) = LA + AL -2e® e. 
(This relation is readily verified; we do not need the dyadic in the sequel.) 



;i.33) 



One need only inspect (1.33) to see that A and L play symmetric roles. In the one- 



charge case it is possible to diagonalize the self-adjoint matrix L in such a way that A 
takes essentially the same form in terms of suitable variables (which are just the action- 
angle variables). This self-duality is no longer present for N^N^ > 0, however. Indeed, 



A is then still self-adjoint, whereas L is not. Even so, (1.33) can again be used to great 



advantage in the construction of the action-angle map, as detailed in [31]. 

We are now prepared to specialize the above to the arena with which the present paper 
is concerned. First of all, we choose 

N+ = 2N, = N. (1.34) 

This choice ensures that there exists a 2A^-dimensional Poincare-invariant submanifold 



1.35) 



with V given by (1.15), of the 6iV-dimensional phase space {^l,uj). This submanifold 
will be described in detail in Section 4, and we shall presently add a brief qualitative 
description. 

Consider now the solution 

exp{tH - yP)Q, Q = {x,p) G fi, (1.36) 

to the joint Hamilton equations for H and P. Such a joint solution exists and is given 



(1.36) as 



by (|1.36|), since the H-fiow and P-flow commute. Using (1.2) and (1.16), we can rewrite 

exp(Mo(MS'_ —vS+))Q = {x{u,v),p{u,v)). 

Mo = V3, c = 7r/3, 



Next, specializing to 



;i.37) 



1.38 
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we introduce 

Tn{u,v) = det{l3N + exp{in{l - 2n)/3)A{x{u,v))). (1.39) 
For a given point Q in the phase space Q, this yields a function depending on n G Z and 

{u,v) e M^. 

We are finally in the position to state the principal result of this paper: When the tau- 



function (1.39) is restricted to the 2A^-dimensional Poincare-invariant subspace flp of fl, it 



coincides with the Tzitzeica tau-function ( 1.14 ) evaluated on P ( 1.15 ). Our demonstration 



of this equality uses in particular special cases of the fusion identities obtained in [33]. 
(For completeness we add the proof of these specializations in Appendix A.) The reduction 
in the matrix size from 3A^ x 3A^ to 2N x 2N hinges on all points in f2p yielding pairs of 
complex-conjugate eigenvalues exp(^j±z7r/3) for the Lax matrix L on the boundary of the 
allowed angular sector. For each of these pairs there is an additional eigenvalue exp(6'j), 
so that the spectrum of L on Qp involves only N degrees of freedom 6i, . . . ,6^ & K. (and 
not the 3A^ of the full phase space). These variables may be viewed as action variables, 
and there are N canonically conjugate 'angle' variables gi, . . . , gA? G M. As a consequence, 



Qp can be identified with V, cf. (1.35) and (1.15) 



A better understanding of how Qp arises as a submanifold of the 6A^-dimensional phase 
space can only be achieved by invoking a great many details concerning the action-angle 
map constructed in [31], which are summarized in Section 4. At this point we only add a 
few qualitative remarks, so as to render these details more accessible. 

First, the above coordinates are not quite action- angle coordinates. Rather, the precise 
definition ofQp involves the harmonic oscillator transform. This transform is an extension 
of the action-angle transform, which takes into account that Q contains an open dense 
subset that is the union of submanifolds for which the angles vary over M^^^' x T', where 
T ~ (— TT, tt] and I takes all values in {0, 1, ... , A^}. From a physical point of view, these 
submanifolds can be regarded as the subsets of Q on which I particle-antiparticle bound 
states are present. Now when the bound state internal actions converge to their minima, 
the /-torus collapses to lower-dimensional tori in precisely the same way as for a harmonic 
oscillator Hamiltonian Yl]=iip'j + ^j)) "which motivates the terminology. 

In these terms, Qp amounts to the subset that arises from the submanifold with N 
bound states by taking the torus to a point, so that each of the pairs is in its ground 
state; moreover, each of the remaining particles has action and angle variables that are 
paired with those of the bound states, in such a way that the asymptotic space-time 
dependence of \E' is that of 3-body clusters moving apart, each of the clusters staying 
together as in the A^ = 1 case. 

In Section 5 we present a close-up of the latter case. For A^ = 1, various questions 
can be rather easily answered, and this special case is also useful as an illustration of 
the general case. In particular, we obtain some explicit information on the 2-dimensional 
space Qp for c G (0,7r/2), and study the Tzitzeica 1-soliton solution corresponding to 
c = 7r/3. 

In Section 6 we compare the Tzitzeica solitons under consideration to the ones obtained 
by Kaptsov and Shanko [16] and to the solitons arising by the above two-step reduction 



( 1.11 )-( 1.12) from a class of 2D Toda solitons constructed via Darboux transformations. 
At face value, these two types of Tzitzeica solitons seem different from the ones obtained 
from the Kyoto solitons in Section 2. As we show, however, the latter are a subclass of 
the former. 
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We have added Section 6 primarily because it yields an affirmative answer to the two 
natural equality questions at issue, but as a bonus it yields a new insight on the Tzitzeica 



tau- function tq (given both by (1.39) and (1.14)): It is in fact positive, and equal to the 
square of a simpler tau-function occurring in the Kaptsov-Shanko work [16]. 

In Section 7 we consider further aspects of our results, in the form of several remarks. 
Specifically, we show that the space-time translation generators restricted to Qp give rise 
to a reduced integrable system, we introduce and comment on space-time trajectories for 
the Tzitzeica solitons, isolate the difference with the setup of [33], and comment on an 
eventual quantum analog of the correspondence between 3-body clusters and solitons. 

A final remark concerns the solitons of the Demoulin system of equations [39,40]. 
As it turns out, their tau-function form is obtained by taking c = it/ 6 in Sections 2 
and 3. It is a challenging question whether they can also be tied in with the relativistic 
Calogero-Moser systems. 

We have relegated a proof of the fusion identities used in Section 4 to Appendix A. 
In Appendix B we collect some auxiliary results concerning pfaffians, which we need in 
Section 6. Finally, Appendix C contains a sketch of the Darboux type construction of 



explicit tau-function solutions to the 2D Toda equation (1.10). 



2 The soliton reduction 2D Toda Tzitzeica 

Our starting point is the 2D Toda M-soliton solutions introduced by the Kyoto school [36, 
37] in the form 

E n e^p(/^.4n)- n (2-1) 

Ati,...,MM=0,l l<j<M l<j<k<M 



jk 



[aj - ak)ibj 



j,fc = l,...,M, 



(2.2) 



{aj - bk){bj - akY 

oo 

e,,. = + nln(a,/6,) + i Y^ikM - ^5') + ^'.-(«7' " ^7'))' (2-3) 
The sequence of 'times' ti^± encodes the hierarchy of commuting flows. Restricting atten- 



tion to the 2D Toda equation (1.10), one should set 

ti,+ = ti^- =0, / > 1, ti,+ = 



(2.4) 



to obtain a solution. 

We proceed in several steps to show how a suitable specialization of the 3M parameters 



a, b and gives rise to a tau-function with the '-Boo '-symmetry ( l-H )• This also involves 
the choice 

M = 2N, N e N*, (2.5) 
which will be in force from now on. 



Proposition 2.1. Define a Cauchy type matrix C with elements 



Cjk = ^^-^, J,A; = 1,...,M, 
ttj - bk 



(2.6) 



8 



and a diagonal matrix 

Dn = diag(exp(^i,„), . . . , exp(^M,n))- 
Then r„ may be rewritten as 



(2.7) 

(2.8) 



Proof. Using Cauchy's identity (1.24), the principal minor expansion of the rhs of (2.8) 
yields 



M 

1=0 /c{i,...,Af},|/|=« jei j,kei,j<k 



which amounts to (2.1). 



(2.9) 



□ 



Next, we specialize the parameters so as to achieve the Boo reduction. Specifically, we 
choose 



-aj, j = l,...,M, 



exp(6) = exp{(pj)aM-j+iiaj + au^j+i) , j = l,...,N. 



(2.10) 
(2.11) 



a^viili-j+i) = - exp{ipj)aj{aj + aM-j+i) \ j = 1, . . . , iV. (2.12) 



Proposition 2.2. With (2.10)-(2.12) in effect, Tn has the B^o symmetry (1.11). 
Proof. Setting 



with 



Xj = exp{{ipj + ipj)/2), j = 1, . . . , M, 

V^Af-j+i = V2j, j = l,...,N, 
i)j = -iv{aj + aM-j+i) + iu{aj^ + a^j_j^^), j = 1,...,M, 
we clearly have 

XM-j+l=Xj, j = l,---,M. 
Also, the reparametrized matrix -D„C is the product of the matrix 



and the matrix with j/c-element 



J = diag(l7v, -1 



O-AI-j+l 



N 



O'M-j+l J Q-j + O-M-fc+l 

After a similarity transformation r„ may be rewritten as 

Tn = \lM + VAnVJl 

where 



j,fc = l,...,M. 



V = diag(xi, ...,Xn,Xn,---, Xi) 

(-a,)"(«M-fc+i)-"+' 



ttj + au-k+i 



j,k = l,...,M. 



(2.13) 

(2.14) 
(2.15) 

(2.16) 
(2.17) 

(2.18) 

(2.19) 

(2.20) 
(2.21) 
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From this definition of An we readily infer that 

A.-n+l,jk = ~A.n,M-k+l,M-j+l- 

Hence we obtain 



(2.22) 



(2.23) 



where the superscript denotes the transpose, and TZm denotes the M x M reversal per- 
mutation matrix. Now we have 



TZmT^T^m — T^i TImJTI'M — — J ■ 
From these formulae we deduce 

r_n+i = I 1m - V'RMA\nMVJ\ = I Im + -RmVA^VJ-Rm \ 
= \1m + VAiVJ\ = |1m + JVAnV\ 



(2.24) 



Tr,. 



as required. 



(2.25) 

□ 



Finally, we show how by a further specialization of parameters, periodic reductions 
of the solitons of the B^o Toda lattice arise. We trade the 2N parameters oi, . . . ,aM 
for N parameters 6i, . . . ,6n and another parameter c, which corresponds to the coupling 
parameter in (1.16)-(1.19): 



ttj = exp{9j — ic), j = 1, . . . , N, 
aM^j+i = - exp(6'j +ic), j = l,...,N. 



Proposition 2.3. With (2.26)-(2.27) in force, we have an implication 

C = Vr// ^ Tn+l =Tn, I = 1,2,.... 



(2.26) 
(2.27) 

(2.28) 



In particular, the tau-function is 3-periodic for c = vr/S, and hence is a tau-function form 
of the Tzitzeica N-soliton solutions. 



Proof. Using also (2.11)-(2.15), we obtain 



exp(e.,n) = ;T^X?e-'^"^ 
2t sm c 



exp(^ 



M-j+l,n) 



2i sin c 



2 2mc 



X^j = exp{(pj - 2sin(c)[t;e''^ + ue ^^]), 



where j = 1, . . . , N. From this (2.28) is obvious. 



(2.29) 

(2.30) 

(2.31) 

□ 
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3 Explicit form of the tau-functions 

In this section we first make the sohton formulae described in the last section more 
explicit, exemplifying them for the simplest case. A suitable reparametrization of the 
(pj^s in (2.11)-(2.12) then yields a sum formula for the arbitrary-c tau-function that is 



not only compact and informative, but which can also be tied in with the relativistic 



Calogero-Moser tau-function (1.39) for the special c-value 7r/3 



First, from (2.1)-(2.4), the 2-soliton solution of the 2D Toda lattice is 

(ai - 62) (oi - as) 



(3.1: 



After the reduction to the B^o Toda lattice described in Proposition 2.2| (|2.3|)-(|2.4|) yield 

exp(0,n) = (-^)X^Af-. + l ^^,K.)^ exp(eM-,+l,n) = 

where 



J - M-j+l o,.(u,v) 
^ 1 



Vj{u, v) = ipj + i[{a- ^ + aj^l_-_^-^)u - {aj + aM-j+i)v] 



(3.2) 
(3.3) 



and j = 1, . . . , N, and (2.2) becomes 

{dj — ak){aM-j+i — ttM-k+i) 



jk 



{aj + aM-k+i)iO'M-j+i + Ofc) ' 



(3.4) 



where j,k = 1,...,M. The 2-soliton solution of 2D Toda (3.1) becomes the 1-soliton 
solution of the B^o Toda equation, given by 



r, 



. _ T , (-l)"(a>2 X) ,yiK^,) , QiQ2(ai -a2)^ ^2^,(„,^) 



fai + a2j 



(3.5) 



The result of the specialization (2.26)-(2.27) is that (3.2)-(3.3) become 



exp(^. 



3(1— 2n)ic 



,ipjiu,v) 



exp(^ 



2t sm c 2t sm c 

iPj{u, v) = ipj -2 sin(c)[ne~^^' + fe^^]. 



-{l-2n)ic 



,iPj(u,v) 



(3.6) 
(3.7) 



where j = 1,...,N. Using (2.26)-(2.27) in (3.4), we finally obtain the following key 
expressions for the functions fjk- 

smh\{d, - 9u)/2) _ cosh^((g,-gfc)/2)-sin^c 

Jjk — . , 2///, n \ . 9 5 Jj,M-k+l — , 2///1 n \ ' 



sinh^((% -4)/2) + sin^c 



cosh^((% - 4)/2) 



for j, = 1, . . . , and 



fjk — fM-j+l,M-k+l, (3.9) 

for j,k = 1, . . . , M. Furthermore, the specialization of the 1-soliton solution of B^o Toda 



(3.5) gives 



= 1 + 



^ ^ cos((2n-l)c) ^^^(,^,) _ cos^c 



zsmc 



4 sin^ c 



(3.10) 
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So far in this section we have simply restated various formulae in the last section in 
more explicit form. Now we make a final reparametrization of the phase constants (pj. 
Specifically, choosing from now on (g, 6*) G V (cf. (1.15)), we set 



exp(v?j) = 2i sin(c) exp(gj/2)Fj(6'), j = 1, . . . 

k=l,...,M,k^j 



(3.11) 
(3.12) 



Substituting this in (3.6)-(3.7), we obtain 

exp(^j-„) = exp((l - 2n)ic + qj{u,v)/2)Fj{6), 

exp{^M~j+i,n) = exp((2n - l)ic + qj{u, v)/2)Fj{e), 
qj{u, v) = Qj - 4sin(c)[fe^^ + ue~^^], 
where j = 1, . . . , N. Note from (3.8) that the functions fjk are positive and hence so are 



(3.13) 

(3.14) 
(3.15) 



Fi, . . . , F/v. 

We have now arrived at the final form of the tau-function. In the proof of the next 
proposition we obtain an explicit sum formula to show it is real-valued. 



Proposition 3.1. The tau-function Tn given by (2.9) with the factors defined by (3.8) 



(3.9) and (3.12)-(3.15) is real-valued. In particular, for c = n/S it yields a real-valued 



N-soliton solution of the Tzitzeica equation. 



Proof. It remains to prove that r„ is real. It is clear from (2.9) that r„ is equal to a sum 



of terms Ts associated with all subsets 5* of {1, ... , M}. To obtain an explicit expression 
for Tg, we write the subsets as 



where 



5= {«!,..., i„M-jb + l,..., M-ji + 1}, 



l<ii<--- <is<N, l<ji<---<jb<N, s,be{0,l,...,N}. 



(3.16) 



(3.17) 



Then we deduce from (2.9) that the contribution to r„ corresponding to 5* is given by 

Ts = exp (^{b - s){2n - l)zc + ^ gi,(n, t;)/2 + J2q,^{u,v)/2^ Fs{e), (3.18) 



(7 = 1 



where 



Fs(0) 



n tf'^- 



Next, defining the reversal permutation r of {1,2, 



follows from (3.9) that 



(3.19) 

.M} by r(j) = M-j + 1, it 

(3.20) 

(3.21) 

where the over bar denotes the complex conjugate. Hence, if r{S) is equal to S then Ts is 
real, and otherwise Ts + Tr(s) is real. Thus r„ can be written as a sum of real terms. □ 



and then from (3.18) that 



T 



r{S) 



Ts 
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To illustrate this result, we write the 1-soliton tau- function in the explicit form de- 
scribed in the above proof. We have 

T% = 1, r{i,2} = exp(gi(M,t;))F{i,2}, (3.22) 

T|i| = exp((l-2n)«c + giKt;)/2)F|i}, (3.23) 

T|2} = exp {{2n - l)ic + qi{u, v)/2) F|2}, (3.24) 



where (recall (3.8) and c G (0,7r/2)) 

F{2} = F{,y = f^^"^ = 1/cosc, F|i,2} = 1. (3.25) 
Clearly, the terms and T{i^2} are real, and so is the sum 

_ _ 2 cos((2n — l)c) / , X , X / X 

T|i} + T|2} = ^ exp gi M, /2 . 3.26 

cose 



The 1-soliton tau-function is given by 

^ ^ ^ 2 cos((2n-l)c) ^^^ ^ 27) 



cose 

In particular, for the Tzitzeica case c = 7r/3 we have 

To = n = 1 + 2e^i("''')/2 + e^^("''') = (1 + e'?i("'^)/2)2^ ^3^28) 

and 

r2 = 1 - 4e«i("''')/2 + e^iC"'^). (3.29) 

Finally, we sketch how to obtain the 2-soliton tau-function in a similar way. In this 
case M equals 4 and for the 16 subsets of {1, 2, 3, 4} there are 10 real- valued combinations 

^0, 7^{l,4}, T|2,3}, 7^{1,2,3,4}, T|i}+T|4}, T|2} + T|3}, 

^{1,2} + ^{3,4}, ^{1,3} + ^{2,4}, ^{1,2,3} + ^{2,3,4}' ^{1,2,4} + ^{1,3,4}, 

whose sum gives the 2-soliton tau-function. 

4 Tzitzeica solitons vs. relativistic Calogero-Moser 
dynamics 

In the Introduction we have already sketched in general terms how the 2A^-dimensional 
Poincare-invariant subset f2p of the phase space uj) arises. In order to fill in the details 
of this qualitative picture, we begin by specifying the subspace of VL on which the maximal 
number N of bound states is present, choosing the internal actions at first larger than 
their minimum, so that there are N angles varying over T^. Moreover, until further notice 
we work with c G (0, 7r/2), since this eases the notation and adds insight on the Tzitzeica 
case c = 7r/3. With this starting point understood, the variables g, 6* in the equations 
(4.20)-(4.24) of [31] should be specialized as follows. 



First, as already mentioned (cf. (1.34)), the particle and antiparticle numbers A'^-i. and 
of [31] must be chosen equal to 2N and (so that the number N used in [31] becomes 
3A^), and the parameters of [31] should be specialized as 

IX = 13 = I, T = g/2 = c. (4.1) 
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Second, the bound state number / of [31] should be taken equal to A^. Thus the numbers 
fc+ and k_ in (4.20)-(4.24) are equal to N and 0, resp., and we obtain 2N particle variables 



gi,...,g;v gM, ^^iv<---<^i, (4.2) 
and 4N bound state variables 

• • • 5 92Ar, ^'at+i, . . . , 6'2Ar e C, ^{On+j) & {0, c), Onj^j 9N+k, (4.3) 

q2N+j = Qn+j + y (1 - i-f), = On+j, (4.4) 

where j, k = 1, . . . , N and k ^ j. 

One advantage of these variables is that they give rise to compact expressions for the 
symmetric functions of the dual Lax matrix A. Specifically, they are given by 



1^1 

with 



5K^) = $^exp(5^g,)p,(^), 

\i\=i \iei J 



l = 0,...,3N, (4.5) 



M(»= n l-''^<'.''^-''^)/^)+-'-l'^^ (4.6) 

Sinh((6'inin(jj) - 6'inax(j,i))/2) 



cf. (5.68) in [31] with parameters (4.1). 



The action-angle variables are now given by 

Xj = qj, p- = e,, (4.7) 



X 



3 



Ij = ^(Qn+j) + {N + j- l)n, 6, = -2^(9^+^), (4.9) 

where j = 1,. . . ,N and 7j G (— 7r,7r] (mod 2tt), cf. (4.26) in [31]. To handle the limits 
Sj I —2c, for which the 7j-torus reduces to a point, we should switch from 7j,5j to the 
'harmonic oscillator variables' Uj,Vj, the minimum 6j = —2c corresponding to the origin 
Uj = Vj = 0, cf. Chapter 5 in [31]. We shall presently analyze the behavior of Si{A) 
under this limit. Before embarking on this rather technical issue, however, we clarify how 
the limit gives rise to a 2A^-dimensional submanifold that is invariant under the Poincare 
(inhomogeneous Lorentz) group. 

To this end we use the explicit description of the commuting flows generated by Hamil- 
tonians of the form 

Hh = T:rh{\nL), (4.10) 

with h{z) an arbitrary entire function (cf. (6.79) in [31]), which follows from (6.6) in [31]. 
It reads 

exp(^if;,)(x^p^; 0, 0) = (xj + th'{pl), . . . , + th'{p%),f; Xi + m{h\p^/2 - zc)), 

. . . , xjv + m{h\pN/2 - ic)),p, 0, 0). (4.11) 
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The time and space translation generators H and P (cf. (1.16)), for which h{z) equals 
Mocosh(z) and Mosinh(z), resp., reduce to 

N N 

H = Mo ^(cosh(pp + 2 cos(c) cosh(pj/2)), P = Mq ^(sinh(pp + 2 cos(c) sinh(pj/2)), 

(4.12) 

on the submanifold at issue, yielding flows 

e*^/^°(x',p'';x,p,0,0) = {x\ + t sinh(p^), . . . , + t sinh(p^),/; Xi +tcos(c) sinh(pi/2), 

. . . , xtv + t cos(c) sinh(pAr/2),p, 0, 0), (4.13) 
g-yP/Moj^^s^^s.^^^^g^g^ = {x\ -ycosh.{p\), . . . ,x% - ?/cosh(p^),/; Xi -ycos{c) cosh(pi/2). 



,xn - ycos(c) cosh(pjv/2),p, 0,0). 



(4.14) 



The boost generator is given by (cf. (1.16)) 



N 



i? = 5^(x^^ + 2x,). 



(4.15) 



whence we have 

e"-^(x^p^;x,p,0,0) = (x^Pl + a, . . . ,p% + a; x,pi + 2a, . . . ,pn + 2a, 0,0). (4.16) 



From these formulae it is obvious that the 2A^-dimensional submanifold Qp obtained 
by requiring 

Xj = cos{c)Xj, pj = 2pj, j = 1,...,N, 
is Poincare invariant. On Qp we have 



(4.17) 



AT 



N 



N 



H = p{c)Mo^coshp'j, P = p{c)Mo^smhp'j, B =p{c)^x'j, (4.18) 



where the prefactor p{c) is given by 

p(c) = 1 + 2 cos(c) 
Also, the symplectic form u reduces to 



N 



UJ 



p{c) ^c?x^ A dp^j. 



(4.19) 



(4.20) 



Thus, when we reparametrize VLp with variables 

qj = p{c)Xj, Oj =pj, j = 1,...,N, 
then q, 9 are canonical coordinates. Moreover the space-time flow reads 



(4.21) 



e^"~y^{q,e) = {qi+p{c)Mo{tsmhei-ycoshei),.. . , g^+j9(c)Mo(t sinh^^-y cosh ^^v), ^), 

(4.22) 
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or, equivalently, 
where 



AT. 



(4.23) 



(4.24) 



With (1.38) in force, this coincides with the Tzitzeica sohton space-time dependence 
(3.15). Of course, this state of affairs is still far from a proof that the tau-function defined 
at the end of Section 3 and the tau-function (1.39) eval uate d on are equal. 

In order to demonstrate this equality, we return to (4.5)-(4.6) and the variables q^ Q. 
On f2p the latter specialize as 



^3 = Qj/Pi^ 



N+j) 



COS 



c)qj/pic 



J 



Oj = ej, 6 



N+j 



e, + tc, j = i,...,iv. 



and (4.4) still holds. More precisely, we have 

G2N+j 



9j — ic, 



1 iV, 



but we can only infer 



QN+j + q2N+j = 2cos{c)qj/p{c) + —(I 



\N\ 



(4.25) 
(4.26) 

(4.27) 
(4.28) 



since the imaginary parts are undetermined in the collapsing torus limit. 

We proceed to analyze this limit for (4.5)-(4.6). First, whenever / contains the index 
N + j, but not the index 2N + j, with j = 1, . . . , A^, or vice versa, then it is clear from 
(4.6) that the limit of pi{9) vanishes. Thus we need only consider I that either contain 
both indices or neither. As a consequence, we only encounter q^+j and q2N+j in the 
combination (4.28), which confirms that the limit is well defined. 

Introducing the breather sets 

B, = {N + j,2N + j], j = l,...,Ar, 
we proceed to simplify the limit of pi{9Y ^oi I of the form 



(4.29) 



(4.30) 



with (|3.17|) in effect. This boils down to a quite special case of the fusion procedure 

(4.31) 



detailed on pp. 237-238 of [33]. Thus we put 

J= {zi,...,^„ji + A^,...,j6 + Ar} c {l,...,2Ar}, 



and 

to obtain 



\2 



limp/(^^)^ 



Oj, Cj = c, Cj+N = 2c, j = 1, . . . , A^, 



j-j- sinh2((r/^- - r]k)/2) + sm^{{cj + Ck)/2) 



sinh ((r^j- - rik)/2) + sin {{cj - Cfc)/2) 



(4.32) 



(4.33) 



cf. (2.27) in [33]. To render this step in our reasoning self-contained, we present a proof 
of (4.33) in Appendix A. 
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Requiring (1.38) from now on, the following two cases arise for the product factor pj^ 
on the rhs of (4.33): 



j,k < N or j, k > N, with j ^ k ^ pjk 
j < N,k > N or j > N,k < N ^ Pjk = 



_ smh\{9j - 9k) /2) + 3/4 
sinh2((^^, -^fc)/2) ' 

sinh2((0j- -0fc)/2) + 1 
sinh2((^^^.-^,)/2) + l/4' 



Taking positive square roots of the pjk, we deduce from (4.33) that we have 



1/2 



limp/(^) = sj n Pjk 

jeJ,kfJ 



(4.34) 
(4.35) 

(4.36) 



for a certain sign sj. We shall determine this sign later on, cf. (4.45) 



Next, we focus on the principal minor expansion of the tau-function (1.39), restricted 



to the submanifold of with bound states present, with internal actions Si, . . . ,Sn in 
(-27r/3,0). It reads 



3N 



Tn{u,v) = y^exp(i7r/(l - 2n)/3)Si{A{x{u,v))) 



1=0 
3N 



^exp(i7r/(l - 2n)/3)^exp{^qi{u,v))pi{e). 

\i\=i 



(4.37) 



(=0 \i\=i iei 

Taking the actions to their minima — 27r/3, we deduce from the above that the limiting 



tau-function is a sum of nonzero contributions Cj for all J of the form (4.31), with 

s 

Cj = exp{i7T{s + 2b){l-2n)/3)Y[exp{qi^{u,v)/2) 



xHexp g,,(«,.)/2 + ^(l-(-)^) n P^i^- 

0=1 ^ ^ jeJMJ 



(4.38) 



We are now prepared to state and prove the principal result of this paper. 



Theorem 4.1. Let c = ir/S. Then the relativistic Calogero-Moser tau-function (1.39) 



restricted to the subspace fip c:^ V is equal to the Tzitzeica tau-function (2.9) with the 



factors given by (3.8)-(3.9) and (3.12)-(3.15) and with parameters in V . 

Proof. We need only show equality of the general contribution Cj to the general term Ts 



in (2.9), cf. (3.18). First, we compare fjk, given by (3.8)-(3.9) with c = tt/3, to pjk given 
by (4.34)-(4.35). From this we easily deduce 



1/2 



Fs{e)= n p% 



(4.39) 



Next, we note that the factors involving qj{u,v) are in agreement. Hence the asserted 
equality comes down to an equality of the remaining numerical factors. Specifically, it 
remains to show 

(-)^xpf^(l-(-)^)Vi = l. (4.40) 
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To this end we now calculate the sign sj i n (4.36 ). We begin by noting that for 
9 E given by (4.2 the product pi{9) (4.6) has a positive numerator. (Indeed, 



each radicand is either positive or has a nonzero imaginary part; in the latter case it is 
matched by a factor with the complex-conjugate radicand.) We are therefore reduced to 
analyzing the phase of 



n 



Sij, 



Sij 



sinh((^jnin(jj) ^max{ij) 



)/2), 



(4.41) 



for the special 6 under consideration, namely, 
9 = {9i, . .. ,9n, 9i+iTc/3, . . . , 9^+^/3, ^i-ztt/S, . . . , 9^-^/3), 



< 9i. (4.42) 



We already know from (4.36) that this phase is just the sign sj. Indeed, we have 



sinh((^j -9k- #)/2) sinh((0j - 4 + #)/2) > 0, ^ = vr/S, 27r/3, 
which confirms that 11/ is either positive or negative. 



(4.43) 



We continue to analyze the contributions to sj from indices in / (4.30), taking (4.42) 



into account. First, we consider the contribution of v ^ I- Due to the ordering of the 



9j in (4.42), any j ^ I with j < N yields Si^j > 0. Also, any that is not a subset of 
/ yields a contribution of the form (4.43) with ip = 7t/3. Hence all indices ii, . . . ,is in I 
yield positive signs. 

Next, we study the contribution of Bj^ C /. For i E I with i < N it again follows 



from (4.43) that we get a positive sign. Now suppose B^ is not a subset of / and consider 



the pertinent product 



^jl3+N,k+NSjp+N,k+MSjij+M,k+NSjp+M,k+M- 



(4.44) 



Both for k < jp and for k > jp this product is easily seen to be negative. Therefore any 
pair Bj , Bj. with Bj included in I and B^ not included in I gives rise to a minus sign in 
11/. Now / contains b breather index sets, so there are N — b breather sets not contained 
in /. Thus we finally deduce 

SJ = (-)^(^-^). (4.45) 



With this explicit formula in hand, it is routine to verify (4.40). Consequently, our 
proof of the asserted tau-function equality is now complete. □ 



5 The N = 1 case 

In this section we consider various aspects of the special case = 1, which yields an 
illuminating illustration of the above constructions and equalities. More specifically, we 
focus on features of the Tzitzeica 1-soliton solution and the space f2p for c G (0, 7r/2). 
For iV = 1 it follows from ([3^-([3^ that we have 



with 



To = ri = 1 + 2F + F2, T2 = 1-AF + F^, F = exp{q{u,v)/2), (5.1) 
q{u,v) = q - 2V3{ue~^ + ve'^) = q + 2V3{t sinh{9) -ycosh{9)), (5.2) 
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cf. (1.2). Thus the two r2-zeros for F = 2± ^/S yield two parallel space-time lines 



y^{t) = ttanh(^) + (q - 2 ln(2 ± v^)), (5.3) 

2v3 cosh(6') 

where \l/ diverges, cf. Fig. 1. In this figure and in Fig. 2, we plot for clarity rather 
than and truncate — at a finite cutoff value. The regions between the singularities, 
where T2 is negative so that \l/ = ln(r2/ro) takes complex values, are approximated by the 
'plateaux' in these figures. 




Figure 1: One-soliton solution. — \E'(t,?/) is shown, the y-axis points left to right and the 
t-axis points into the page. The 'plateau' approximates the region in which ^^(t, is 
complex- valued. 



Recalling (1.39), we see that the r2-zeros correspond to x{{t, y) or x^it, y) being zero. 
The space-time lines may therefore be viewed as the trajectories of the two particles in 
the 1-soliton cluster. 

For = 1 it might seem an easy matter to locate the 2-dimensional space Vtp within 
the 6-dimensional phase space VL. In fact, however, it is not even easy to find the 1- 
dimensional submanifold oiVtp consisting of if-equilibrium points, i.e.. 



(5.4) 



in explicit form. (Note that for c = 7r/3 this yields the stationary 1-soliton tau-functions. 



Of course, it is immediate from ( 1.16 )-( 1.19) that we need 

Pi =P2 =Pi = 0, 



(5.5) 



for H to have an equilibrium. Requiring this, one expects from physical considerations 
to get an equilibrium point for 



0, 



(5.6) 



and a suitable d. We now confirm this for c G (0,7r/2). Then we have from (1.16)-(1.19) 



i/(c/, -rf, 0, 0, 0, 0) 



2(/+r 



vl/2 



+ r, r 



1 + 



sin^ c 



sinh^ d ' 



r 



1 - 



sin^ c 



cosh^(rf/2) 



• (5.7) 



19 



Thus we have H 3Mq for c? ^ cxd and iJ ^ oo for c? ^ 0. Since H is equal to 
Mo cosh(^j), it has an absolute minimum Mqp{c) for 6 = {0,ic, —ic). Now we can 
still choose xl,xi G M (cf. (4.7)-(4.8)), so H has a 2-parameter family of stable equilibria. 
It is straightforward to check that for d equal to 



c?e = cosh ^(1 + cosc), 



we have 



1 + 2 cos c 
(2 + cose) cose' 



r 



[1 + 2 cose) cose 

2 + cos e 



Hence the rhs of (5.7) equals Mop(e), so we do get a stable equilibrium for d 
implies that all of the points 

E(or) = (4 + cr, -4 + cr,cr,0,0,0), a G M, 



(5.^ 



(5.9) 



4- This 



(5.10) 



are also equilibria. 

Since H equals Mop(e) on Qe, one might expect that the 1-parameter family of equi- 
libria -^(0") yields Qe- In fact, however, only E{0) belongs to Qe- To be specific, we 
assert that it corresponds to g = 0. To show this, we use (4.5) with = 1 to find the 
symmetric functions of A on Qe- We determined the relevant limits below (4.28). In 
particular, (4.33) yields p(e)^ for / equal to {1} or {2, 3}, since = on Qe- Using also 



(4.25) and (4.28), we readily obtain 



Si{A) = p{c)exp{q/p{c)), 

S2{A) = -p(e)exp(2cos(e)g/p(e)), 
53(A) = -exp(g). 
For g = this implies that the spectrum of A is given by 

a{A) = {l + cos(e) ± (cos^(e) + 2 cos(e))^/^ -1}. 



(5.11) 

(5.12) 
(5.13) 



(5.14) 



The first two eigenvalues can be written as exp{±de), cf. (5.8). Thus the origin of Q.e 



corresponds to generalized positions xj" = de,X2 = —de,Xi = 0. Since ^Ie consists of 
if-equilibria, we also have pt = pt = Pi = 0. Hence we deduce i^^(O) & Qe and g = 0, as 
asserted. 

On the other hand, the equilibrium points inQp with g 7^ are harder to find explicitly. 
As already announced, they do not include the equilibria E{a) for a 7^ 0. (Indeed, E{0) 
does belong to Qe, as just shown. Now the cx-shift of the generalized positions corresponds 
to a shift of xf and Xi by a, so that the f2p-condition xi = cos(e)x^ no longer holds true for 
cr 7^ 0.) Rather, they are obtained by acting with the commuting flow exp{yP/Mo), y G M, 
on the equilibrium E{0), yielding translated equilibria 



T{y) = (x+ (y) , x+ (y) , xr (y) , 0, 0, 0) , ye 



In view of (5.13) and (4.22) this yields the sum rule 

xfiy) + ^tiy) + xiiy) = p{c)y. 



(5.15) 



(5.16) 
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Also, from Hamilton's equations for P/Mq we have 

a;+(yy > cos(c), j = l,2, x'^ {y)' > cos\c) , (5.17) 

so that Xi{y),X2{y) and x^{y) are strictly increasing functions of y. 

Next, we deduce from (5.11 )-(5.13) that we have the reflection symmetry 

xt{-y) = -^tiy)^ Xii-y) = -Xiiy), (5.i8) 

so it remains to determine the functions for y > 0. This is presumably possible, but 
we have not pursued this. However, the large-y asymptotics can be established from the 
above. Specifically, setting 

e = exp(cos c — l)y, 



we obtain for e — > 



xt{y) = y + lnp{c) + 0{e^), 
xtiy) = cos{c)y - - Inp(c) - -^/p{c){l - p(c)"^)e + O(e^), 

x~{y) = cos{c)y- ^lnp(c) + ^^/p(cj{l - p{c)~^)e + 0{e^). 



(5.19) 

(5.20) 
(5.21) 

(5.22) 



A final observation of interest concerns the kinetic and potential energy density of the 
stationary 1-soliton solution, i.e., the functions 

EKiy) = idy^)V2, Ep(|/) = exp(^) + exp(-2^)/2 - 3/2, ^ = ln(r2/ro), (5.23) 



obtained from (5.1) for = 0. It is far from obvious, but true that these functions 
are equal. This equality can be verified directly by a straightforward, but quite tedious 
calculation. 

A more conceptual derivation of this virial type identity will now be given. First, we 
note that any t-in dependent solution to (1.1) satisfies the ODE 

-2/ 



yy 



— e 



Now from (5.1) we obtain a solution to (5.24) of the form 
f{y) = (7(exp((g - 2v^t/)/2)), g{z) = In 

Defining 
it satisfies 



1 + 2^ + ^2 



2± v^. 



e5(") > 0, 



Z > Z+, Z < Z- 



and 



exp(^(l)) = 



■ 0, z ^ ±oo, 
-1/2, ^7'(1)=0. 



(5.24) 

(5.25) 

(5.26) 

(5.27) 
(5.28) 

(5.29) 
(5.30) 
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Next, consider the Hamiltonians 



H± (x, p) = p72 - 14 (x) , 14 (x) = ±e^ + e-'72 - 3/2 



(5.31) 



The potential —V^{x) has a maximum at x = and yields a Newton equation 

X = - (5.32) 



Comparing this to ( |5.24[ ), (5.27) and (5.29), we see that g{z) for z > Zj^ corresponds to the 
E = {] orbit with x{t) < coming from for t — > — cxo, and g{z) for z < z^ to the E = {) 
orbit with x(t) < going to for t oo. By energy conservation, we have x^/2 = 14 (x), 
which implies Exiy) = Ep(y) for the y-intervals where C'^^^-* > 0. 

It remains to show the identity for the y-interval where e^^^^ < 0. Then we can 



compare (5.24) to the Newton equation 



X 



-e - e 



-2x 



(5.33) 



corresponding to The E = orbit stays to the left of the origin and has its turning 
point at X = ln(l/2), cf. (5.31 ). Comparing this to ( |5.30 ), we see that it corresponds to the 
real part of g{z) for z G (z_,z+). Hence EK{y) = Ep{y) now follows from x^/2 = \4(x). 



6 The Darboux and Kaptsov-Shanko solitons 

In this section we begin by showing that the solitons obtained by the i?oo-reduction 
from the Kyoto 2D Toda solitons are equal to those obtained by a similar reduction 
from a seemingly different class of 2D Toda solitons. The latter are constructed via 
repeated Darboux transformations. As will transpire, this procedure yields a larger class 
of solutions, involving an arbitrary constant antisymmetric matric C. In order to obtain 
equality to the solitons of Section 2, this matrix must be suitably specialized. 

Our demonstration of equality leads to a new representation of the latter solitons. 
This representation can be exploited to show that tq equals the square of a simpler tau- 
function. This is because it readily leads to tq being the determinant of an antisymmetric 
M X M matrix A. Hence it follows that we have 



^0 



Pf(A). 



(6.1) 



The pfaffian can be explicitly evaluated, and when the Tzitzeica substitutions of Sec- 
tions 2-3 are made in A, then the resulting r is the one obtained by Kaptsov and Shanko 
in their study of the Tzitzeica equation [16] . 



We proceed with the details. We start from the tau-function r„ (2.19), with J', V and 



A given by (2.17), (2.20) and (2.21), and with the quantities Xj ^ arbitrary at this 
stage. We claim that r„ is equal to the determinant fn of the matrix 



where 



^^-l-7^^ ) 

B = diag(/3i, . . .,Pm) 



(6.2) 

(6.3) 
(6.4) 
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■njk 



dj + ak 



n+l 



j, = 1, 



,M, 



provided that j3i, . . . , Pm are chosen such that 



AT. 



(6.5) 



(6.6) 



(Recall TZn denotes the N x N reversal permutation matrix.) 

In order to prove this, we denote the columns of Ain by Ci, . . . , cm, and use 

Tn = \Mn\ = |Col(Ci,. .. ,Cm) I = | Co1(Cm, • • • , C7V+1 , -c^v , . . . , -Ci) | . (6.7) 

Thus, f„ equals the determinant of the matrix 



-aj)"'{aM-k 



1m + J^nJ' , J^n,jk — f^jf^M-k+r 

ttj + ttM-k+l 

Transforming this matrix with the similarity matrix 

S = diag(7i, . . . , 7a/), 7i = iPM-j+i/Pj)^^'^, j = 1 

we obtain 



-n+l 



,...,M, 



with 



fn = \lM + VAnVJl 

V = diag((A/3M)'/^ {(32Pm^iY'\ {PmPiY^'). 



Hence the /3-constraint (6.6) entails equality of r„ and t„, as advertised. 
We can now compare the new representation 

Tn = \Cn + VAnV\, 

obtained by choosing 



Xj, j = l,...,N, 



(6.8) 

(6.9) 

(6.10) 
(6.11) 

(6.12) 
(6.13) 



(so that (6.6) is obeyed, and B and V both equal V), with the solitons obtained by a -Boo 
symmetry reduction from the Darboux type 2D Toda solitons. (See Appendix C for a 
sketch of their construction.) The reduced solitons are of the form 

T^ = \C + BAnB\, (6.14) 

where C is an arbitrary antisymmetric M x M matrix and B is a diagonal matrix 

^ = diag(/3i,...,/3M), (6.15) 

with diagonal elements of the form 



Pj = aj exp{—ivaj + iua- '^), j = 1, 



M. 



(6.16) 



Recalling the definition (2.13)-(2.15) of Xj, we see that (6.6) is obeyed when we choose 
ctM such that 

(6.17) 



ajaM-j+i = exp(v9j), j = 1,...,N. 



23 



Provided we also specialise the arbitrary antisymmetric matrix C to C-ji (cf. (6.3)), we 
therefore conclude equality of the reduced Darboux type solitons to the reduced Kyoto 
solitons. 

We continue by using the new representation (6.12) to show the square property of tq, 

(6.18) 



cf. (6.1). As it stands, the matrix C-ji + VAqD is not antisymmetric. But we have 

1 



0,jk 



ttj + ak 2 



1 



a,- — Ofc 



aj + ak 



so that we can write 



A = ^c®c-^A c = (i, 



where A is the antisymmetric matrix with elements 



A.^ 



■jk 



ttj — ttk 

Ctj -\- Ctk 



As shown in Appendix B, the determinant 



To = \Cn + -{VC(»VC-VAV)\ 



equals the square of 



r = Ff{Cn-VAV/2), 
and the pfaffian can be explicitly evaluated. The resulting formula is 

N I 



T 



E E n<'i.. n 



c 



Jmjn ' 



with 



1=0 l<ji<-<jt<N m=l 
1 



l<m<n<l 



, . . . , X T , 



Cij — AijAi^M-j+lAj^M-i+lAM-j+lM-i+i^ h 3 — ■ ■ ■ y 



cf. (B.21). 



Next, we specialize ai,. . . ,aM as in ( 2.26 )-( 2.27). Using (2.31), this yields 

cos(c) 



di 



2i sin(c) 



exp{(pj — 2 sin(c) [vkj + uk- ]) 



2 ,.2 



+ Kj +2 COs{2c)KiKj 
Hi + l^j J Kf + H?, —2 COs{2c)K,iKj 



where 



Kj = exp( 



J = l, 



AT. 



(6.19) 



(6.20) 



(6.21) 
(6.22) 

(6.23) 

(6.24) 
(6.25) 



(6.26) 
(6.27) 
(6.28) 



We now substitute (3.11), and then choose c = ir/S and introduce new parameters 

kj = VSkj, exp{sj) = exp{qj/2)Fj{e)/2, j = 1,...,N, (6.29) 

-V, y = —u. (6.30) 



X 
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Then we finally obtain 



dj = exp(sj + xkj + Sykj "^), 

h-kj V kf + k'j -kjkj 
ki ~\~ kj J kf -\- kj -\- kikj 



(6.31) 
(6.32) 



With (g, 6) varying over V (1.15), we have kj > and exp(sj) > 0, so that dj > and 



Cij > 0. The function r is therefore positive, implying tq is positive. Note that r can be 



rewritten in the form (2.1), with M, ^j^n and fjk replaced by A^, dj and Cjk, respectively. 



Last but not least, the tau-function (6.23) with the substitutions (6.31)-(6.32) co- 



incides with the tau-function obtained in Section 2 of [16]. Moreover, combining the 
relations 

exp(^) = Ts/ro, To = Ti, tq = r^, (6.33) 



(cf. (1.13), (1.11) and (6.1)), and the 2D Toda equation of motion (1.10) with n = 1, we 

(6.34) 



deduce 

2dudJnT = 1 -exp(^). 
Therefore, the function exp(\l/) coincides with the function v employed in [16]. 



7 Concluding remarks 

(i) {Integrability on Qp) It is not hard to see that the flows generated by the Hamiltonians 

Hr = Tr(L"), r G M*, (7.1) 

leave Qp invariant, provided 



cos rc 



cos c . 



(7.2) 



Indeed, this readily follows from (4.11 ) and (4.17), noting that (7.1 ) corresponds to (4.10) 



with h{z) = e^^. Thus the restricted phase space {flp, J2f=i d'lj ^ d9j) and the Hamilto- 



nians 



Hr, ±r = 1 + 27rk/c, k e Z, 



(7.3) 



give rise to an integrable system on Qp. For the Tzitzeica case c = 7r/3, one can choose 
as the N independent Hamiltonians the power traces 



Tr(L 



= 0,l,...,iV- 1. 



(7.4) 



(ii) {Space-time trajectories) As we have shown in Theorem 4.1, on Qp the tau-function 



(1.39) equals the Tzitzeica tau-function, and hence is real, cf. Proposition 3.1. This reality 



property is far from obvious for n = 0, 1, but for n = 2 reality on all of Q is in fact clear 



from real-valuedness of A, cf. (|1.32|). Specializing to Qp, we have 

2A'" 

T2{U,V) = 



N 



JJ[1 -exp(a;+(M,f))] ■ JJ[1 + exp(Xj. {u,v))]. 



(7.5) 



1=1 



Now it follows from Section 6 that to{u,v) is positive. Therefore, the Tzitzeica solution 

^{u, v) = Ht2{u, v)/Toiu, v)), (7.6) 
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has logarithmic singularities at the zeros of T2, as we have already seen for = 1 in the 
previous section. We proceed to analyze these in terms of the space-time coordinates t 



and y, cf. (1.2). 



Fixing t, the function T2 is positive for y ±00, and has 2N sign changes for finite y. 
The locations of these zeros on the y-axis are distinct for all t (since < ■ ■ ■ < x^). 
Thus one obtains 2N space-time trajectories, which may be viewed as the locations of 
the particles in the A^-soliton solution. For t ±00 these trajectories exhibit soliton 
scattering, with a factorized phase shift in terms of the function In(cjj), cf. (6.32). A more 



systematic analysis would be feasible by following the path laid out in Chapter 7 of [31], 
but this is beyond our present scope. See, however. Fig. 2 for a plot of the 2-soliton 
collision. 




Figure 2: Two-soliton interaction. — \E'(t,y) is shown, the y-axis points left to right and 
the t-axis points into the page. The 'plateau' approximates the region in which '^{t,y) is 
complex- valued. 



(iii) {Comparison with [33]) The 2D Toda soliton tau-functions studied in [33] do not 
include the above real-valued Tzitzeica tau-functions. This is because in [33] the condition 
r_„ = r„ is imposed (cf. (2.10) in [33]), whereas for the Tzitzeica case we have r_i = T2 ^ 
Ti = Ti. To accommodate this different starting point, the fusion procedure in Section 2 
of [33] starts from tau-functions that in terms of the relativistic Calogero-Moser systems 
amount to 



Tr, 



\1m + exp{-2inc)A{x{u,v))\, A{x) = diag(exp(a;+), . . . , exp(a;|^)). 



(7.7) 



Even so, we could specialize the fusion identities of [33], since they only pertain to the 
action variables, and the dependence on the latter is governed by the same function ((2.15) 
in [33]) for particles and antiparticles. 

(iv) {Quantum analogs) We have chosen A^_|_ = 2A^ and = N throughout the paper 



(cf. (1.34)), but we could just as well have started from particles and 2A^ antiparticles. 



Indeed, this amounts to working with the 'charge conjugate' dual Lax matrix 

A^{x) = diag(exp(x|'), . . . , exp(x^), - exp(xi"), . . . , - exp(x2^)), (7.^ 
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and tau-f unction 



(m, v) = det(l37v - exp(z7r(l - 2n)/3)A^{x{u, v))), (7.9) 



instead of (1.32) and (1.39). Clearly, one can then proceed in the same way as before. 

On the other hand, one cannot enlarge the 3-body-soliton correspondence without 
venturing into complexified phase spaces, losing control of the action-angle map in the 
process. More precisely, there exist real-valued Tzitzeica soliton tau-functions that corre- 
spond to two different charges, hence giving rise to breather-like bound states. (Indeed, 
this is not hard to see from the explicit form of the solitons at the end of Section 3; for 
example, one need only perform a suitable analytic continuation in the 2-soliton solution 
to obtain the 1-breather solution.) However, it can be shown that these do not correspond 
to subspaces of the real Calogero-Moser phase spaces with arbitrary and N^. 

It may be expected that this picture persists on the quantum level. To be specific, a 
suitable reduction of the unitary joint eigenfunction transform for the commuting quantum 
Hamiltonians with c = tt/3 and 3A^ variables (which, to be sure, is not even known to 
exist to date) should give rise to a unitary transform with variables, which can be 
interpreted as a transform for N quantum solitons with the same charge. However, no 
such reductions are likely to exist when different charges are involved. In particular, if 
quantum breathers for the Tzitzeica quantum field theory do exist (a feature that is taken 
for granted in most of the work within the form factor program) , then they are not likely 
to have analogs in the Calogero-Moser particle picture. (By contrast, for the sine-Gordon 
case a complete correspondence is expected [41].) 

(v) {Demoulin solitons) From Propositions 2.2 and 2.3 it is clear that for 

c = 7r/6 (7.10) 



the tau-function satisfies 
As a consequence, we have 



ri = To, T2 = T5, T3 = n. (7.12) 



Setting 

hn = Tn+lTn-l/rl = exp(0„ - 0„-l), (7.13) 



(where we used (1.9)), we deduce 

hi = ho, h2 = /i5, h^ = /i4, (7.14) 



and 



Now from (1.3) we have 



Hence, setting 



hih2h3 = 1. (7.15) 
{\nhn)uv = '2hn-hn+i-hn-i. (7.16) 

h = hi = T2/Ti, k = h3 = T2/T3, (7.17) 
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and using (7.15), we obtain 



(In/i), 



h 



1 

hk' 



(Ink), 



k 



1 

kh' 



(7.18) 



This system of relativistic wave equations is the Demouhn system, cf. [40], p. 343, 
Eq. (9.58). Therefore the c = 7r/6 tau-functions with parameters in V yield real- valued 
Demoulin solitons in tau-function form. In particular, from (3.15)-(3.27) we see that the 
one-soliton case is given by (7.12) and 

rl = l + 2G + G^ r2 = l + G^ rg = 1 - 26* + G^ G = exp{q/2 - ve^ - ue'^). (7.19) 



Appendices 



A Fusion 



In this appendix we detail how the formula (4.6) for pi{9Y leads to (4.33) in the collapsing 
torus limit. We begin by recalling that in this limit pi{9) vanishes, unless the subset / of 
{1, . . . , 3A^} has the form 

I={jh (A.l) 

where 

J = Jn{l,...,2iV}, /, = {j}, In+j = {N + j,2N + j], j = l,...,iV, (A.2) 



cf. the paragraph preceding (4.30). For a subset / of this form, (4.6) entails 

p0f 



n n 



sm\i\{dm-k)/'2)+sm^c 



sinh^((^„-^„)/2) 



(A.3) 



Next, we substitute (4.26) and (4.27) in this formula and cancel terms in the interior 



product using fusion identities. A general fusion identity, in which the sets Ik are of 
arbitrary cardinality, can be found in [33], but here we only detail the special cases we 
need. 

Using the trigonometric/hyperbolic identity 

sinh^ X + sin^ y = sinh(x + iy) sinh(a: — iy), (A. 4) 



the general pair factor in (A.3) can be written as 



sinh((g^ - g„)/2 + ic) sinh((g^ - g„)/2 - ic) 
sinh2((^^ - e^)/2) 



(A.5) 



For a soliton-soliton interaction included in J, there is only one pair factor in the product. 



Specifically, we have Ij = {j}, Ik = {k} and (from (4.26) and (4.32)) 6j = rjj, Ok = rjk, so 
the contribution to the product for a soliton-soliton pair is 

sinh^((?7j - r]k)/2) + sin^ c _ smh{{r]j - rik)/2 + ic) smh{{r]j - r]k)/2 - ic) _ s{2)s{-2) 



sinh ((?7j - r^fc)/2) 



sinh ((?7j - rik)/2) 



s2(0) ' 
(A.6) 
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where we have set s{l) = smh.{{rij — r]k + lic)/2). For a sohton-breather interaction we 
have Ij = {j}, I 
yielding two pair factors 



J - us, ^N+k = {N + k,2N + k} and 9j = r]j, e^+k = Vk + ic, 92N+k = Vk - ic, 



s(l)s(-3) s(3)s(-l) _ s(3)s(-3) _ smh^{{rij - rjk) /2) + ?>m^{'ic/2) 



s{l)s{ 



sinh^((r7j - r/fc)/2 + sin2(c/2) 



(A.7) 



Finally, for a breather-breather interaction in J, we get iN+j = {N + j, 2N + j}, iN+k 
{N + k, 2N + k}, and hence four pairs whose contribution to the product is 

s(2)s(-2) s(4)s(0) s(2)s(-2) s(0)s(-4) s(4)s(-4) smh^{{r]j - %)/2) + sin2(2c) 



s2(0) s2(2) s2(0) 



sH-2] 



s2(0) 



sinh {{r]j -r]k)/2) 



(A.8) 

Using rji+N = Vi Q = c, c^+i = 2c for / = 1, . . . , iV (cf. (4.32)), these three cases may 
be encoded in the single formula 



sinh^(fa- - r/fc)/2) + sin^((9 + Cfc)/2) 
sinh^ {{r]j -r]k)/2 + sin^((9 - Cfc)/2) 



j, k = l, 



,2N. 



(A.9) 



Hence (4.33) results. 



B PfafRan identities 



In this appendix we show that the determinant on the rhs of (6.21) equals the square of 
the pfaffian on the rhs of (6.22), and that the latter has the explicit form (6.23)-(6.25). 
We proceed in a slightly more general way, since this eases the notation and adds insight. 

First we note that the determinant is of the form |f + A|, where A is a 2A^ x 2A^ 
antisymmetric matrix. Now we invoke the following lemma. 

Lemma B.l. Assume A is an antisymmetric 2N x 2A^ matrix and v G C^^. Then one 
has 

\v®v + A\ = \A\. (B.l) 
Proof. By continuity it suffices to prove this for an invertible A. Then we have 



\v®v + A\ = \A\\l2N + A'^v ® v\ = \A\{1 + {v,A"^v)). 



(B.2) 



Now since {A 



-i\t 



A ^, it follows that A ^ is also antisymmetric. Hence the 

□ 



inner product {v,A ^v) = v*'A vanishes, and (B.l) results. 
As a consequence, we have 

To = \Cn-VAV/2\=T^, 



(B.3) 



with r given by (|6.22|). We continue to show (|6.23|). Again, we first study a slightly more 

(B.4) 



general situation. We begin by recalling the definition 



-A, 
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where the sum is over all index choices with 

1 = < ■ ■ ■ < <2N - 1, ii<ji,...,iN< Jn, 
and a is the permutation 

a : {iuji, . . . ,^jv, Jat) ^ (1,2, . . . ,2A^). 
Denoting Pf (v4) by [1,2,..., 2A^], this implies an expansion formula 



(B.5) 



(B.6) 



2N 



[l,2,...,2iV]=^(-)"[l,n][2,...,n,...,2iV], 



(B.7) 



n=2 



where the notation on the rhs will be clear from context. Consider now the pfaffian of 
the antisymmetric matrix 



B 



A7^ 



N 



-A7^7v 
The relevant (upper triangular) part of B is 
Au ■ ■ 

An-i,n An-i,n+i a + An-i 



+ A, AgC. 



(B.8) 



N+2 



X + A 



N,N+1 



A 



N,N+2 



A 



N+l,N+2 



2N 



(B.9) 



^2Af-2,2Af 
A2N-l,2N 



Therefore we have 



N 



Pf(5) =^A^-'5;, 



(B.IO) 



/=0 



where Si is a sum over certain 'minors'. Specifically, taking (B.7) into account to check 
signs, we readily get 



Sk= Yl [-^'i' • • • ^Jk^^N-jk + 1, . . . ,2iV - ji + 1]. 

l<ji<-<jk<N 



(B.ii: 



(Hence in particular 5*0 = 1, S2N = [1, • • • , 2A^].; 



Applying this expansion to r (6.22), we obtain 

N 



/ 1 \ 

^ = E E n (-2^1 ■^■i,---,Jz,2iV-j, + l,...,2iV-ji + lU, (B.12) 

1=0 l<ji<---<ji<N m=l ^ ^ 

where [...]^ denotes the reduced pfaffians for A. We now obtain an explicit formula for 
the latter. 
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Lemma B.2. Suppose A G M2i(C) is given by 

ctj — Clk 



A 



■jk 



aj + ak ' 



j,k = l,...,2l. 



Then we have 



i<i<fe<2i 



(B.13) 



(B.14) 



This explicit evaluation is proved in [42]. It is presumably known, but we do not know 
a reference. Indeed, up to its sign this lemma can be viewed as a consequence of Cauchy's 
identity. To see this, note that the lemma entails 



1^1= n ^ 



2 

■jk- 



(B.15) 



^<j<k<2l 

If we now substitute aj exp(2gj), then ( B.15[ ) becomes 

|(tanh(gj - qk))\2ix2i = JJ tanh^(gj - q^), 

i<i<fc<2« 



(B.16) 



an identity that readily follows from Cauchy's identity (cf. (B42) in [43]). Thus, Cauchy's 
identity implies (B.15), so (B.14) follows save for its sign. 

Returning to (B.12), we substitute ( |B.14[ ) in the ^-pfaffians, yielding 



n 

j<k,j,k&{ji,...,ji,2N~ji+l,...,2N-ji+l} 



Introducing 



we can write (B.17) as 



(B.17) 



Cjk — AjkAj^2N-k+lAk,2N-j+lA2N-k+l,2N-j+l, ^ j < k < N, (B.18) 



Substituting this in (B.12) and setting 

1 



d] = -^XjA,2N-j+i, j = 1,...,N, 



we finally obtain 



N 



I 



^=5Z J2 n 

1=0 l<ji<-<ji<N m=l j<k,j,ke{ji,...,ji} 



which amounts to (6.23). 



(B.19) 



(B.20) 



(B.21) 
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C Solutions obtained by binary Darboux transfor- 
mation 

The Aoo Toda lattice (1.3) has a family of solutions obtained by an iterated binary Dar- 
boux transformation. We follow the derivation given in [44], but make some small mod- 
ifications so as to conform with the notation and conventions used in the rest of this 
paper. 



A Lax pair for (1.3) is given by 



ic.i] 



(Thus, equality of cross derivatives entails (1.3).) It is covariant with respect to the 
Darboux transformation [45] 



4>n = (pn + ln(a„/a„_i). 



(C.2) 
(C.3) 



where a„ is an eigenfunction (i.e., any particular solution of (C.l)). The adjoint Lax pair 
is 



- Bn,u = i exp(0„ - 0„_i)5„_i, Bn,, 



iB. 



n+1 



+ 



(C.4) 



Given an eigenfunction a„ and an adjoint eigenfunction bn for a seed solution 0„ = 
0° = ln(r°_,_]^/r°) to (1.3), an eigenfunction potential f2(a„, bn) is defined, consistently and 
uniquely up to an additive constant, by the three requirements 



«exp(0°+i - 4>l)bnan 



Vt(cLn-, bn\u ~ 

Q{an,bn) - Q{an,bn-i) = -b^a^. 



+1) 



(C.5) 

(C.6) 
(C.7) 



The Lax pair (C.l) is also covariant with respect to the binary Darboux transformation 

an^{An,bn)/^{an,bn), 



A„ — Ar, 



(Pn = (Pn + ln(^^(a„, 6„)/il(a„_i, &„-i)). 



(C.8) 
(C.9) 



or, equivalently, = f2(an-i, &n-i)T°. (In (C.8) An denotes the general solution to (C.l) 
with replaced by 0°.) 

The Mth iterated binary Darboux transformation is expressed in terms of an M- 
vector An, whose components satisfy (C.l) with 0°, and an M- vector B„, whose 
components satisfy (C.4) with cfp. These eigenfunction and adjoint eigenfunction 
vectors are used to define an M x M matrix eigenfunction potential VLn = fi(A„,B„,) 
satisfying 



^n,u = -iexp(0„_^i 



A 



n+l, 



-iBn 



+1 



An,, fiji ^n—1 



The solutions obtained are then expressed as 



det(fin-i)i"^ 



Bn ® An. 

(C.IO) 
(C.ll) 
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From now on we omit the " and write f„ as r„. 

The Boo Toda lattice has 0_„ = — 0„ for G Z, so that (po = and the system 
is semi-infinite, expressed in terms of 0i,02,---- This reduction can be achieved by 
choosing solutions satisfying r^n+i = Tn- In this reduction, it is consistent to choose 
adjoint eigenfunctions given by B,„ = (— 1)"A_„. It then follows from (CIO) that 

= (C.12) 

and so provided M is even, M = 2N say, we indeed obtain r^n+i = Tn- 

In particular, to obtain soliton solutions for the A^o Toda lattice, one chooses as seed 
solution the vacuum solution = 1,0^ = 0, and then vectors of eigenfunctions and 
adjoint eigenfunctions given by 

A„,,- = a,-a-"e-™'^^+*"'^7\ B„,, = j = 1, . . . , M, (C.13) 

where aj,(3j,aj and bj are arbitrary constants. Then, one obtains 

Cik + Pi^k ' ^ , e-^"(°''-^^)+^"("fe''~''^"'^ , (C.14) 
' ' ak-b, J 

where Cjk are arbitrary constants. In the i?oo -reduction, one must take 6j = — and 
Pi = at giving 

Tn = det ( C,k + a^.a, t^^l^^^Le-™(%+-^)+-K-'+«7')^ . (C.15) 
\ ttj + ttk J 

Here, Cjk are arbitrary constants satisfying, in accordance with the symmetry condition 
( |C12[ ), Ckj = —Cjk- 

Thus we have arrived at the expression (6.14) that is the starting point of comparison 
with the Kyoto solitons in Section 6. 
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